RECONSTRUCTION OF THE TRANSMISSION COEFFICIENT 
FOR STEPLIKE FINITE-GAP BACKGROUNDS 



IRYNA EGOROVA AND GERALD TESCHL 



Abstract. We consider scattering theory for one-dimensional Jacobi opera- 
tors with respect to steplike quasi-periodic finite-gap backgrounds and show 
how the transmission coefiicient can be reconstructed from minimal scattering 
data. This generalizes the Poisson-Jensen formula for the classical constant 
background case. 



1. Introduction 

In classical one-dimensional scattering theory the transmission coefficient can 
be reconstructed from the reflection coeflicient via the well-known Poisson-Jensen 
formula. This formula plays a crucial role in inverse scattering theory since it shows 
how to compute the left scattering data from the right one and vice versa. Moreover, 
it is also one of the key ingredients for deriving the associated sum rules which 
have attracted an enormous amount of interest recently (see e.g. [17], [21], [23] . 
[24j . [31]). Furthermore, these sum rules are intimately connected with conserved 
quantities of the associated completely integrable lattices (see [2^, [22 )• Finally, 
the reconstruction formula can be viewed as the solution of a scalar Riemann- 
Hilbert factorization problem which arises in the nonlinear steepest decent method 
[5] when deriving the long-time asymptotics (see [13], [18], respectively [T9]). 

Moreover, the same is true in case of scattering with respect to a finite-gap 
background [6], [30]. In this situation the analogous formula was given in [28] 
including the associated sum rules (see also [7], [S], [H]). Again there is a close 
relation with the solution of a scalar Riemann-Hilbert factorization problem on the 
underlying Riemann surface which arises in the nonlinear steepest decent method 
[H, [H], [16], and [20]. 

However, while scattering theory with a steplike constant background is classical 
and goes back to the early sixties [3] (see [5] for the most recent results), even in this 
case the reconstruction formula is unknown to the best of our knowledge except for 
the case when the two spectra overlap. This might be related to the fact that the 
case where the two spectra do not overlap will not be solved in terms of elementary 
methods, but will already require tools from the theory of elliptic surfaces, as we will 
see below. In case of steplike finite-gap backgrounds, scattering theory is again well- 
understood by now [5], ^U\, however, the reconstruction formula is again unknown 
to the best of our knowledge (except for the case of two finite-gap backgrounds in 
the same isospectral class [8]). The main purpose of our paper is to fill this gap 
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and provide a reconstruction formula for the left/right transmission coefficient in 
terms of the left /right scattering data. 

2. Notation 

We begin by introducing some required background from the theory of hyper- 
elliptic curves to be used in the remainder of this article. For further information 
and proofs we refer for instance to [1] , [3] , [11] , [12] , or [27] . 

Let M be the Riemann surface associated with the function P^/^(z), where 

23+1 

(2.1) P{z)= lliz-E,), Eo < El <■■■ < E2g+u 

3=0 

g E N. M is a compact, hyperelliptic Riemann surface of genus g. We will choose 
P^/^(z) as the fixed branch 

25+1 

(2.2) Pi/2(z) = - n V^^' 

j=o 

where ^ is the standard root with branch cut along (— oo,0). 

A point on M is denoted hy p = {z, ±P^/^{z)) = {z, ±), z G C. The two points 
at infinity are denoted by p = oo±. We use 7r(p) = z for the projection onto the 
extended complex plane C U {oo}. The points {{Ej,0),0 < j < 2g + 1} C M are 
called branch points and the sets 

g 

(2.3) n± = {(z,±pi/2(z)) I zeC\I]}cM, j:=\J[E2j,E2j+i], 

3=0 

are called upper and lower sheet, respectively. Note that the boundary of II-i- 
consists of two copies of E corresponding to the two limits from the upper and 
lower half plane. 

Let {ajjbjj^j^i be loops on the Riemann surface M representing the canonical 
generators of the fundamental group 7ri(M). We require aj to surround the points 
i?2j"-i, E2j (thereby changing sheets twice) and bj to surround Eq, £'2j-i counter- 
clockwise on the upper sheet, with pairwise intersection indices given by 

(2.4) Gj o Uk ^ bj o bk ^ 0, ajobk = Sjk, ^<j,k<g. 

The corresponding canonical basis {Cj}^=i for the space of holomorphic differentials 
can be constructed by 

(«) £ = £*)^, 

where the constants c(.) are given by 
The differentials fulfill 

(2.6) / Cfc = <5j,fc, / Ck = Tj^k, 7-j,fe=rfej, l<j,k<g. 

J aj J bj 

For further information we refer to [TT], [?fl App. A]. 
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In addition, we will need Green's function (in the potential theoretic sense) of 
the upper sheet 11+ : 

Lemma 2.1 ([28|). The Green function o/H-)- with pole at zq is given by 
(2.7) g{z,zo)^-Ref uip^p^, p = (z, +), po = (zo, +), 

JEo 

where po — po* (i.e., the complex conjugate on the other sheet) and ujpq is the 
normalized Abelian differential of the third kind with poles at p and q. 

Clearly, we can extend g{z, zq) to a holomorphic function on M\{po} by dropping 
the real part. By abuse of notation we will denote this function by g{p,po) as well. 
However, note that g{p,PQ) will be multivalued with jumps in the imaginary part 
across ^-cycles. We will choose the path of integration in C\[£^0: ^2g+i] to guarantee 
a single- valued function. 

From the Green function we obtain the Blaschke factor (cf. [28] ) 



(2.8) Bip, p) = exp (g{p, p)) = exp ( Upp^j, tt{p) e R, 

which has the following properties: 
Lemma 2.2. The Blaschke factor satisfies 

(2.9) B{Eo,p) = l, and B{p\ p) ^ B{p, p*) ^ B{p, p)-^- 

it is real-valued for tt{p) G (— oo, Eq). 
Moreover, 

(2.10) \B{p,p)\ = l, pes, arg(B(p,p)) =5,(p), 7r(p) G [iJa.-i, ^2,], 
where we set E^i = — oo, i?2g+2 = oo and 

fO, J = 0, 

(2.11) 5M=\\k^Pf>'^ J = 

[o, J =.9 + 1. 

Proof. The first part including the fact that \B{p,p)\ = 1, p G E, is proven in [28l 
Lem. 3.3]. To see the formula for the argument first observe that Wpp* is real- valued 
on 7r^^(R\I]) and purely imaginary on 7r^^(S]). This can be seen from the explicit 
expression (|2.13p for Upp* given below. Hence, taking the path of integration along 
the lift of the real axis, we see that the integral is real for p = (A,±) with A < i?o 
or A > i?2g+i. For p = (A, ±) with A G [£'2^-1, £-2^] the imaginary part is constant 
and given by half the hj period. □ 

The above Abelian differential is explicitly given by 

/ pl/2+pl/2(p) pl/2+pl/2(^) . rf, 

^ \ 2{n-nip)) 2(7r-.(g)) 

where Ppq{z) is a polynomial of degree g — 1 which has to be determined from the 
normalization uipp* — 0. In particular, 

(2.13) topp, = ( 11!^ + E^^, (^)^ '^^ 



TT - tt{p) pp ' ' ) pi/2 
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3. Reconstructing the transmission coefficient 
Let be two quasi-periodic finite-band Jacobi operators Q 
(3.1) i/±/(n)=a±(n)/(n + l) + a±(n-l)/(n-l) + 6±(n)/(n), fef{Z), 
associated with the hyperelliptic Riemann surface of the square root 

2g± + l 



(3.2) p1/'(z) = - n V^-^f' E^<Et<---<Et,^+i, 

where (7± € N and ^ is the standard root with branch cut along (— oo,0). In 

fact, iJ± are uniquely determined by fixing a Dirichlet divisor J^^j^iP-f^ where 
/i± = inf^af) with fif e [E^j^i,Efj] and af G {-1,1}. The spectra of 
consist oi g± + I bands 

(3.3) a± a(ff±) = |J [E^,,Ef^^,]. 

We are interested in scattering theory for the operator 

(3.4) Hf{n)^a{n-l)f{n~l) + b{n)f{n)+a{n)f{n + l), n G Z, 

whose coefficients are asymptotically close to the coefficients of on the corre- 
sponding half-axes: 

±C30 

(3.5) ~ + - < 

Tl = 

Let ip^{z,n) be the Floquet solutions of the spectral equations 

(3.6) H^ip{n) ^ z^{n), z € C, 

that decay for z e C\a± as n ^ ±cx3. They are uniquely defined by the condition 
tpf{z,0) = 1, ^pf{z,-) e ^^(Z±). The solution 4'^{z,n) (resp. tjj~{z,n)) coincides 
with the upper (resp. lower) branch of the Baker-Akhiezer functions of (resp. 
H-), see [27]. 

The two solutions '>p±{z, n) of the spectral equation 

(3.7) HiP = z^, ze C, 

which are asymptotically close to the Floquet solutions i^^{z, n) of the background 
equations (|3.6p as n — > ±oo, are called Jost solutions. 
Next, we introduce the sets 

(3.8) cr(2) = (7+ n CT_ , cr£^ = cr± \cr(2), (7 = (7+ U Cr_ , 

where a is the (absolutely) continuous spectrum of H and a'^^^ U a'^}^ , cr^^^ are the 
parts which are of multiplicity one, two, respectively. 

In addition to the continuous part, H has a finite number of eigenvalues situated 
in the gaps, ad = {Ai,...,As} C R \ tr (see, e.g., [25]). For every eigenvalue we 
introduce the corresponding norming constants 

(3.9) 7±,i = El'/-'±(^fe'«)l'. l<k<s. 

nez 



^Everywhere in this paper the sub or super index (resp. "— ") refers to the background on 
the right (resp. left) half-axis. 
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Note that this definition has to be slightly modified in the unlikely event that 
ip^{z, n) and hence i^±[z, n) has a pole dX z = \k (see [10] for details). The trans- 
mission and reflection coefficients are defined as usual via the scattering relations 



(3.10) T^{X)ij±{\,n)=^^{\,n) + R^{X)ij^{\,n), A € cr=p. 

Here the values of 'ip±{X, n) for A S ct± are to be understood as limits from above 
'>p±{\ n) = limeio '4>±{^ + i^j n) (the corresponding limits from below just give the 
complex conjugate values '0±(A, n) = lim^jo V'±(-^ ~ i^j "-))• 

The following result is an immediate consequence of pUl Lem. 5.1]. 

Theorem 3.1 ( ,10j ). Suppose a{n), b{n) satisfy (j3.4p . thena{n), b(n) are uniquely 
determined by one of the sets of its "partial" scattering data S-\. or S-, where 



5± = {i?±(A), Aea±; \T±{X)f,Xea^^ 
(3.11) Ai,...,A, eM\((7+Ua_), 7±4>---:7±,. e 



This leads to the natural question if there is a simple way to compute 5+ from 
S- and vice versa (i.e., without solving the inverse scattering problem). It turns 
out that this question reduces to the reconstruction of the transmission coefficient 
T±{z) from S±. In fact, this follows from the following lemma. 

Lemma 3.2 ([[10]). The transmission coefficients T±(z) admit a meromorphic ex- 
tension to C\(T. In general they have simple poles at the eigenvalues Xk of H . In 
addition, there are simple poles at /i^ £ M\a± which are not poles of ip^{z,l) 
(i.e., ~ ^\) and simple zeros at fj,J € R\o'zp which are poles of ify^^z, 1) (i.e., 
aj — ^1). A pole at fi^ could cancel with a zero at fij or could give a second order 
pole if fif ^ Xk. 

Moreover, the entries of the scattering matrix have the following properties: 

(a) p+{z)T+{z) = p^{z)T^{z), 

(b) ^ = i?±(A), Aea«, 

(c) l-|i?±(A)|2 = ^|r±(A)p, Aea(2), 

(d) i?±(A)r±(A)+i?^(A)T±(A) = 0, Aea(2), 

where 

3.12 p±{z) = . 

Hence, the problem is to reconstruct the meromorphic function T_|_(z), z G C\(t 
from its boundary values 

|T+(A)|2, Aea«, 
(3 13) ^|r+(A)p = ^(l-|i?+(A)|2), Aea(2), 

^ = i?+(A), Aea«. 
T+(A) + 

That is, we know its absolute value on (T_ and its argument on the rest cr^^ . There 
will be three Riemann surfaces involved, the one corresponding to a — a+ U a- 
and the ones corresponding to a±. All objects corresponding to a will be denoted 
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as in Section [21 while the objects associated with a± wiU have an additional ± 
sub/supscript. 

Theorem 3.3. The transmission coefficient can he reconstructed from the 

reflection coefficient R+{z) and the eigenvalues Xj via 

VjeM- / \tj.+ eM+ J \k=i J 

(3.14) expf^^/ Qlog(|r+|)c.,,. 

+ / Q(log(^)+log(l-|i?+p))c.... 

27ri J ^(2) ' >+' 

(3.15) + J^^^^ Q ( arg(i?+) + <5-)c.,,-^ , 

where the integrals are taken over the lift, of the indicated spectra to the upper sheet 
IIu (of the Riemann surface associated with a). Moreover, we use the convention 
that we identify z with [z,-\-), and similarly for \k, ^J^f, whenever used in the 
argument of a function defined on a Riemann surface. Here 

(3.16) Af± = {^if \^if eR\c7 and af ^ -1}, 

(3.17) Qiz) = Y\ ^/z — e j , where Cj are defined via \^[e2j,e2j+i] = 

3 3 

and 



with (cf. Lemma \2.2]) 

s 

(3.19) 5i=- ^i(Pj)+ E Siifip + Y^SiiX,). 

Proof. We start by considering the multivalued function 
(3.20) 



t+iz)=i n B^{z,^,J)-'\i n B^iz,^4)\il[B^{z,X,)\T+iz) 

which has neither zeros nor poles on n„ and satisfies 
('|<+(A)|2 = |T+(A)p, Ae^_, 
[arg(t+(A)) = arg(T+(A)) + V, A e n \E-,_^,E-^. 



(3.21) 



Moreover, the absolute value of tj^(z) is single- valued and hence its logarithm is a 
harmonic function on Hi, which can be reconstructed from its boundary values. To 
accommodate the fact that we know its absolute value on (t_ and its argument on 
we consider 

(3.22) Q(z)log(i+(z)). 
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Note that since t+{z) might still have zeros and poles on a, the function log(<+(z)) 
might have logarithmic singularities on a. 

Since Q{X) is real- valued for A € M\f7^^' and purely imaginary for A €E cr^"^ , we in- 
fer that the real part of Q{z) log(t_|_(z)) is harmonic on n„ and can be reconstructed 
from its boundary values 
(3.23) 

Re(Q(A)log(t4A))) = |^(')^r'^"^-(')")' , ' ^ V) 
^ ^ \iQ(A)(arg(T+(A)) + 57), X e n 

using Green's function: 

Re(Q(z)log(i+(z))) Q log(|r+| V,,. 

(3.24) + i ^^^^ Q( arg(r+) + ^ )c.,,. j . 

Dropping the real part we get 



(3.25) exp|^^^^^ Qiog(|T+|)a;,,. -H^^^^^^Q(arg(r+) + 5-)c.,,.^. 

In fact, by [29, Thm. 1] both the left-hand and the right-hand side have the same 
absolute value and hence can only differ by a constant with absolute value one (in 
particular, the right-hand side is single- valued since the left-hand side is). This 
constant must be one since both sides are real-valued for real-valued z to the left 
of CT. □ 
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